Equilibrium spin currents and magnetoelectric effect in magnetic nanostructures 
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We discuss the problem of equilibrium spin currents in ferromagnets with inhomogeneous magneti- 
zation. Using simple microscopic models we explain the physical origin of equilibrium spin currents. 
Next we derive the equilibrium spin current from the Hamiltonian with a gauge field associated 
with local rotations in the spin space. Several examples of magnetic systems are studied in details, 
and the persistent spin current is found to exist in the ground state of these systems. We also 
demonstrate the possibility to measure the equilibrium spin current using the magnetoelectrically 
induced electric field near the ring. 
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1. Introduction. The problem of generation of the pure 
spin currents attracted a lot of attention recently. It is 
mostly related to the perspectives of possible applications 
in spintronics, for which the generation and manipulation 
of spin currents is of a prime importance. 1 Several differ- 
ent ways have been proposed to create the spin currents, 
like, for example, the injection of spin-polarized carriers 
from a magnetic metal or semiconductor^ optical excita- 
tion of spin-polarized electrons or holes j 3 - equilibrium spin 
currents induced by the spiral states in ferromagnets p 
spin Hall effect^ spin transport in presence of the spin- 
orbit (SO) interaction^ and others. 

Recently, the transport of magnetization by magnons 
has been studied by Meier et ali They demonstrated 
that by using a finite-length spin chain between mag- 
netic reservoirs, the pure spin current can be generated 
without the transport of electron charge. Then Schiitz et 
al proposed to use the magnons in a mesoscopic Heisen- 
berg ring under an inhomogeneous magnetic field as a 
possible way to induce the persistent spin current £ It 
was also found in these works^' 8 that the magnetization 
transport by magnons is accompanied by an electric po- 
larization, which can be experimentally observed in the 
vicinity of the magnetic wirei or the mesoscopic ring. 8 
Most recently, Katsura et al used an electronic model of 
two transition atoms mediated by an oxygen atom, to 
show that the equilibrium spin current in this system is 
related to the electric polarization through a new mech- 
anism of the magnetoelectric effect^ which, in fact, is 
much stronger than the previously discussed^ relativis- 
tic mechanism. 

At the same time, the theoretical studies of the mag- 
netization transport provoked a lot of discussions about 
the existence and physical meaning of the non-vanishing 
equilibrium spin currents and about a proper defini- 
tion of the spin current in the non-equilibrium transport 
phenomenaii2*ii*i£ The main idea of these works was to 
re-define the transport spin current, which does not in- 
clude the equilibrium (supercurrent) part responsible for 
the spin torque— 

Here we show that several key points concerning the 
equilibrium spin currents and electric polarization dis- 



cussed in the above-mentioned papersL2i2ii£ should be 
further clarified. In particular, we emphasize that the 
equilibrium spin current (also called supercurrent be- 
cause it is associated with the rigidity of the order 
parameter—) is related to a noncollinear magnetic order- 
ing, and it exists at T = without magnons (in contrast 
to the claims of Refs. 8,12]). The details of electronic 
structure and/or SO interaction are not important for 
the equilibrium spin currents. On the other hand, the 
observable electric polarization is induced by the equilib- 
rium spin current as well as by the magnon transport, 
and the corresponding mechanism necessarily includes 
the spin-orbit interaction. The magnitude of this effect 
strongly depends on the electronic structure of the mate- 
rial, so that it can be much larger than a bare relativistic 
effect. 11 We believe that the effective enhancement of the 
SO coupling for the magnetoelectric mechanism in the 
noncollinear magnets is an analogue of the renormaliza- 
tion (giant enhancement) of the Rashba SO interaction 
in the layered systems, as well as the reinforcement of 
the anomalous Hall effect, which is also related to the 
SO interaction. 

2. Two spins in local fields. We start from a sim- 
ple example demonstrating that the equilibrium spin cur- 
rent does not vanish for two interacting quantum spins 
S = 1/2. Namely, we consider the Hamiltonian with in- 
teracting spins Si and S2 in local magnetic fields Bi and 
B, 



H 



-J Si • S2 — Bi • Si — B2 • S2 



(1) 



where |Bi| = | B2 1 = B. Calculating the time-derivative 
of spin operator Si^ — j- [H, Si^] with (1), we find 

Sifi = — t^vx S\ v S2X + -r £fj,v\ Siy Bi\ . (2) 

The first term in the right side of (2) is the variation 
of spin related to the spin current ii_ > 2i anc ^ the second 
one is related to the local field Bi acting on Si. In the 

equilibrium state, (Si^ = due to compensation of the 

moments transferred by the spin current and generated 
by the field Bi. 
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To calculate the value of spin current in the ground 
state, we use the perturbation theory assuming the in- 
teraction weak, | J | /B <C 1. Let us take the quantization 
axis z along the field Bi, and the axis y in the plane of 
vectors Bi and B2. In the limit of \ J\ /B <C 1, the eigen- 
function of Hamiltonian (1) corresponding to the ground 
state with energy e = —2B is (tp\ = (ipi \ (g> (-i/> 2 | with 



(Vxl = (1, 0) , {fo\ = 



cos-, -1 am- 



is) 



where 8 is the angle between the vectors Bi and B2. 

Calculating the spin current in this state we find 
(jfUa> = -Jsmd/h, and (jf^ 2 ) = (jf^ 2 ) = 0. The sim- 
ilar calculation gives us (jf—i) = — (ii-> 2 )) which means 
nonconservation of the spin current in the presence of the 
local fields Bi and B 2 . The spin current transferred to 
the spin Si generates the torque rotating this spin from 
the orientation along vector Bi. 

3. Hubbard model. This result can be explained by 
microscopic calculation using a two-site electronic model 
with Coulomb interaction. Let us consider a model of 
two-electron system with strong on-site Hubbard interac- 
tion U. This is an example of an electronic system, which 
is often used as a model leading to the localized magnetic 
moments with antiferromagnetic interaction. We show 
that by applying on-site magnetic fields Bi and B 2 with 
different orientations, we induce the spin current trans- 
ferring a spin moment. The Hamiltonian is 



H = -t (c\ a c 2a + c\ a ci Q ) +U ( c iT 



Ci\ Cil 



i=1.2 



(4) 



where t is the hopping, c ia and Ci a are the creation and 
annihilation operators for electrons at the site i with spin 
a. We take |Bi| = |Bi| = B and consider the case of 
strong on-site interaction and weak field, B <C t <C U. 
The operator of the spin at site I is Si = \c\ a cr a pCip^ 

and the spin dynamics is determined by Si = i [H, Si], 
where the dot means time derivative. Using (4) we find 



S f = §(4 



C2P - c\ a c w ) 0% + e^a^c\ a c lf3 . (5) 



Here the first term in the right hand site is the operator 
of spin current transferring spin from site 2 to site 1 . The 
second term is related to the force acting on the spin Si 

by the the field Bi. In the ground state = 0, which 

can be satisfied due to the nonvanishing spin current 



r 



it 



2 CT a/3 \ C la C 2/3 



C L Cl/3 



(6) 



flowing from site 2 to site I. The spin current transfers 
the angular moment compensating the force produced by 
the field B x . 

We consider the system with two electrons, and take 
the basis of two-particle states |TJ.,0), |0,TI), |T>T)> 



Hit}) |l) I}- The Hamiltonian (4), and the operators of 
spin Si and the spin current j can be presented as ma- 
trices in the basis of two-particle states. For example, 
if we take Bi along z, and B 2 in the y — z plane, the 
Hamiltonian in the new basis acquires the following form 
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B + 



H = 



B'J 

(7) 

where we denoted B = B\ z and B'^ — B^^- 

In the limit of strong interaction, t/U < 1, and B = 0, 
this Hamiltonian can be transformed by a unitary trans- 
formation to the block form, in which the part describing 
the spin coupling is H = J(Si-S 2 -l/4)-Bi-Si-B 2 -S 2 , 
where J = 4t 2 /U. By using the same transformation to 
the spin-current matrix, we find that it transforms into 
j = J Si x S 2 . This calculation shows that, microscopi- 
cally, the equilibrium spin current results from the hop- 
ping of electrons between different sites, with an effective 
transfer of moment by the spin of electrons. 

4. Noncollinear ferromagnet. Now we consider a tex- 
tured ferromagnet, which represents a magnetic system 
in a topologically nontrivial metastable state. A simple 
example of such system is a magnetic ring with the easy- 
axis magnetic anisotropy, so that the magnetic moment, 
oriented along the ring, creates a vortex. 

Let us take a continuous classical model of a ferromag- 
net described by the Hamiltonian, which includes the ex- 
change interaction, anisotropy, and the interaction with 
an external magnetic field B(r) (this model was used re- 
cently to determine the Berry phase of magnons in tex- 
tured ferromagnets^) 



H = d A v 



|(^0 2 +^{n(r)}- fiBptin 



(8) 



where n(r) is the unit vector oriented along the magneti- 
zation at the point r, a is the constant of exchange inter- 
action, T {n(r)} is a function determining the anisotropy. 

We will use the definition of spin current related to 
the transformation of the Hamiltonian under local ro- 
tations of vector n in the spin spaced This definition 
is in accordance with a general definition of currents 
in the quantum field theory^ In the case of a classi- 
cal magnetic system with Hamiltonian (8), the corre- 
sponding transformations of the vector n(r) belong to 
the group SO (3). Thus, we perform a local rotation 
n(r) — > R(r) n(r) using the orthogonal transformation 
matrix R(r) = ^W^e^W^e*'^, where V, 0, (f> are 
the Euler angles determining the arbitrary rotations of 
the coordinate frame, and J x , J v , and J z are the gener- 
ators of 3D rotations around x, y and z axes, respectively. 

The Hamiltonian of exchange interaction (the first 
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term in Eq. (8)) in the rotated frame has the form 
H ex = | J d 3 r (din^ - A lpv n v ) 2 , 



(9) 



where the gauge field Aj(r) = (djR)R -1 . The matrix 
Ai(r) can be presented as Ai(r) = iJ^A^ir), where 
Aj(r) belongs to the adjoint representation of the group 
SO (3). Then the exchange energy can be written as 



d, 5, 



a/3 



a/3 ) n & 



(10) 



and the spin current density is defined as j£ = 
7 (5H/8A.i), where 7 is the gyromagnetic ratio. We find 

3i = -iCsJ^np (di S ai - iA\ J s ai ) n 7 , (f f ) 

where c s = 7a. The spin current (11) is gauge invariant. 
We can fix the gauge by taking the auxiliary field A% = 0. 
Then, using the relation «J^ A = e P vx, we finally obtain 



ji = c s n x diii. 



(12) 



Thus, the spin current (12) is nonzero in the noncollinear 
ferromagnets. In particular, it is nonzero in a metastable 
state of the ferromagnet with topological excitations. 16 
This definition is consistent with the spin-conservation 
equation relating the variation of the magnetization in 
time to the divergence of the spin current j&il ri^ + 
divj^ = Tp, where T„ = je^x n v B e x f f , and B e ^ is 
the local effective magnetic field, comprising the local 
external field and other local terms, such as magnetic 
anisotropics. The physical meaning is that of a hydrody- 
namic equation for n(r), where is a local source term 
(torque) . 

The thermodynamic average of the spin current is 



tii) = Cs^iiuX \{n v ) di {nx) + {8n v di 8n\)] 



(13) 



where (X) = Tr (Ie^ fi )/Tr (e~^ H ), and SX = X-{X) 
is the fluctuation of X. Thus, there are two contributions 
to the spin current: the principal one arises from spatial 
variations of the average value of the magnetization axis, 
while the second one is due to the magnetization fluc- 
tuations (magnons). The magnon term is likely to be 
a small correction as compared to the principal (ground 
state) term. The ground state term can be interpreted 
as the supercurrent associated with the rigidity of the 
order parameter Ai&iSi The magnon term has been dis- 
cussed thoroughly^iSiiS and will not be discussed further 
here. In the following, we assume T = 0, so that the 
magnon term vanishes (except for the effect of quantum 
zero-point fluctuations, which are usually negligibly small 
in ferromagnets). In this case, we can drop the angular 
brackets, which means that the thermodynamic average 
is implied. 

5. Mesoscopic ring. To calculate the spin current in 
the ring geometry, it is convenient to use the cylindric co- 
ordinates (p, ip, z) for the point on the ring. We assume 



that the vector n does not depend on p and z. Then the 
exchange Hamiltonian H ex can be written as 



a Co 
2R 



I 



dip ( d v n p + — e^xv $a n v 



(14) 



where R and Co are the radius and the cross section of 
the ring, respectively, &x = LA* t and L = 2nR. 

Let us take Alt constant along the ring. Then $ M is the 
flux of the /^-component of the gauge field Alt through 
the ring. Now we get the azimuthal spin current density 
in the ring as jji — (7/C0) {dH/dQ^). Using (14) and 
taking <5> M = 0, wc find 



■p = ri 



L 



dip {—n z (n p + d v n v ) + n v d v n z ] , (15) 
dip [n z (— n„ + d v n p ) - n p d v n z ] , (16) 



3% = — / dip [n<p {n v ~ d v n p ) + n p (n p + d^n^)} . (17) 



In particular cases of the tangential magnetization 



1, n„ 



0), radial magnetization (n p 



1, 



iitp = n z = 0), and also for any intermediate case with 
the in-plane magnetization making a constant angle with 
the tangent vector, using (15)-(17) we obtain 



2-KCs/L , 



(18) 



It should be emphasized that the spin current (18) is 
related to the assumed metastable state of magnetiza- 
tion field but not to the magnons. The contribution of 
magnons exists for T / but it is small because the 
magnons are weak excitations over the metastable state, 
completely vanishing in the limit of T — > 0. 

The other example is a magnetic ring with uniaxial 
anisotropy in a homogeneous magnetic field along the 
axis z. Due to the anisotropy and exchange interaction, 
the magnetization along the ring is oriented with a cer- 
tain angle out of the ring plane. There exist a metastable 
state a with crown-like magnetization profile in this sys- 
tem. (Due to the dipolar forces it can be the ground 
state.) We can calculate the angle 9 using Hamiltonian 
(8) with ^ 7 {n(r)} = \ n 2 z /2 and A > 0, and the magnetic 
field B oriented along the axis z. Then, using the polar 
coordinates and assuming that n v and n z do not depend 
on the coordinate along the ring, we find the energy of 
the metastable state 



E = tt(oR [anl/(2R 2 ) + \n\j2 - f3Bn z ] 



(19) 



and calculate the angle 6 minimizing the energy (19)^ 
cosfl = max {l, 0B/{\ - a/R 2 )} for A > a/R 2 , and 6 = 
for A < a/R 2 . The spin current in the ground state with 
9 ^ does not vanish. We calculate the components of 
it using Eqs. (15)- (17) 



2-kc s 



2ttc, 



sin ft cost/, j 



sm 2 0, 4=0. (20) 
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Obviously, we obtain the same value of the spin current 
for different magnetization profiles if they can be trans- 
formed to each other by a global rotation. 

6. Electric polarization. As discussed in Refs. 00, 
the spin current implies an electric polarization P, = 
eij^jj /c 2 , which is the relativistic effect of a transforma- 
tion of magnetization to the electric field in the moving 
frame (this fact that has been known long ago22>). The 
important point is that the electric polarization appears 
not only due to the spin transport of magnons but also 
due the spin supercurrent in the ground state (first term 
in Eq. (13)). This was shown by Katsura et al. 9 using a 
microscopic model with two transition atoms and an oxy- 
gen atom in-between. Their model calculations demon- 
strate that the magnetoelectric coefficient determining 
the magnitude of effect depends essentially on the mate- 
rial parameters. Taking it into account, we can write the 
polarization 

P = a me [n (V • n) - (n • V) n] , (21) 

where we denoted by a me the magnetoelectric coefficient. 
This equation implies an enhancement of the bare rela- 
tivistic effect in the condensed matter, i.e., it is due to 
the newly proposed magnetoelectric effect^ It should be 
noted that the suitable materials for the observation of 
polarization are magnetic insulators and semiconductors 
because the induced electric polarization in a good metal 
would be completely screened by free electrons. 

In the cases of the magnetic ring with the tangent mag- 
netization along the ring or the constant radial magne- 
tization, using (21) we find that the polarization vec- 
tor P is oriented along the radius of the ring, and 
P r = —a me /R, where R is the ring radius. 

Let us estimate the magnitude of effect for a typical 
magnetic nanostructure, namely, a ferromagnetic disc of 



radius R and thickness h with a circular vortex domain 
(tangential magnetization) i2i The effect of the vortex can 
be approximated by replacing the disc by a ring of inter- 
nal radius r equal to the radius of the vortex core. In 
this case the polarization P is radial and decays as the 
inverse distance from the ring axis. We can approxi- 
mately estimate the electrostatic potential U at a point 
located on the axis at an altitude z above the disc in the 
limit of r, h <C \z\ -C R. We obtain U ~ — a me h/\z\. 
We estimate a me using the result of a three-atom model 
from Ref. for the polarization per unit volume, which 
gives a me ~ (e/a) (V/A) 3 . Here a is the lattice con- 
stant, V is a hopping energy (hopping between the tran- 
sition and oxygen atoms in y|), and A is a character- 
istic electron energy (between d and p orbitals of dif- 
ferent atoms). This formula was obtained for large SO 
splitting Ago 3> A, and it does not explicitly depend 
the SO interaction. For a magnetic disc, using typical 
values R — 1 /im, ft, = 10 nm, and r = 10 nm, this 
yields U ~ 0.1 (V/A) 3 V at an altitude \z\ = 100 nm. 
The vortex can be removed by application of an exter- 
nal magnetic field. The effect can be possibly detected 
experimentally by using a single-electron transistor as a 
detector of electric field. 

We calculated the equilibrium spin current in meso- 
scopic systems, and showed that it does not vanish in 
the metastable state of magnetic system with topologi- 
cally nontrivial magnetization profile. The electric polar- 
ization induced by the spin current can be experimentally 
measured. 
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